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We investigate the emergence of chiral meso-structures in one-dimensional fluids consisting of 
stacked discotic particles and demonstrate that helical undulations are generated spontaneously from 
internal elastic stresses. The stability of these helical conformations arises from an interplay between 
long-ranged soft repulsions and nanopore confinement which is naturally present in columnar liquid 
crystals. Using a simple mean-field theory based on microscopic considerations we identify generic 
scaling expressions for the typical buckling radius and helical pitch as a function of the density and 
interaction potential of the constituent particles. 
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I. INTRODUCTION 

Chirality is a key signature of life since many biological 
molecules owe their specific function to their molecular 
shape and symmetry. Examples range from molecules 
possessing at least one carbon atom with four different 
substituents to more complicated helical molecular as¬ 
semblies such as the DNA double helix. Apart from the 
microscopic domain, chirality is expressed abundantly 
in the macroscopic world, for example in the left- and 
right-handed helical conformation of snail shells, the spi¬ 
ral patterns of leaves on a plant stem (phyllotaxis), and 
other biological and inert structures. The structure of 
condensed phases are naturally affected by the chiral¬ 
ity of the building blocks. For example, liquid crystals 
composed of chiral constituents may form extraordinary 
structures, ranging from cholesterics, chiral smectics and 
blue phases mm- The particular meso-structure im¬ 
parted by chiral correlations on the micro-scale, such as 
a helical director field in cholesterics or a cubic arrange¬ 
ment of defect lines in blue phases, endow these materi¬ 
als with extraordinary optical and mechanical properties. 
Important applications of chiral mesostructures reside in 
domains such as non-linear optics [3], photonics [4] and 
stereoselective catalysis [5]. 

Helical structures may also arise in systems of achi¬ 
ral particles. Illustrious examples are spontaneous chiral 
symmetric breaking of homopolymers composed of at¬ 
tractive non-chiral monomers [5] , the formation of Bernal 
spiral clusters in patchy colloidal particles [71 [5] and the 
self-assembly of achiral magnetic nanoclusters into helical 
superstructures [9]. Helical meso-structures also emerge 
from packing colloidal particles into cylindrical tubes at 
high pressure nnum where transitions between different 
cluster morphologies may occur upon varying the pack¬ 
ing load [12]. In fact, coiled, toroidal and helical architec¬ 
tures seem to be a recurring feature in systems of (block 
co-)polymers [131 [H] and other soft matter compounds 
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confined in nanopores [THUS]. 

Although the formation of the helical meso-structure 
occurs by means of spontaneous self-assembly, there 
should be no preference for a particular helix sense (left- 
or right-handed) if the constituent molecules are achi¬ 
ral. Spontaneous breaking of the chiral symmetry usually 
happens at the initial stage of the self-assembly process 
by some combination of fluctuations (of e.g. thermal na¬ 
ture) favoring one chiral symmetry over the other. The 
local bias toward a particular chiral symmetry then pro¬ 
liferates to mesoscopic length scales by means of some 
self-sorting and amplification mechanism which secures 
the uniform mesoscale symmetry of the assembly [ninH]. 

In this work we investigate the possibility of helical 
meso-structures in columnar assemblies of discotic par¬ 
ticles with non-chiral pair interactions. Columnar meso- 
phases are characterized by a two-dimensional crystalline 
arrangement of columns, each consisting of a quasi one¬ 
dimensional stack of particles with a disordered inter¬ 
nal structure. In discotic liquid crystals supramolecu- 
lar chirality manifests itself along the columnar direction 
and can be imparted by chiral intramolecular interactions 
m, or by non-chiral ones via some type of directional¬ 
ity in the face-to-face interactions (so-called 7r-stacking), 
or a helical organization of the orientational or transla¬ 
tion degrees of freedom of the discotic particles along the 
backbone of the stack [70] . 

The focus of this paper is to explore the latter type of 
chirality, namely the one expressed by an effective helical 
shape of the columnar stack. To this end we devise a 
microscopic model and use simple mean-field theory to 
identify generic expressions for the free energy associated 
with weak conformational changes from a linear stack to 
a helically buckled assembly of stacked particles subject 
to a confining external field. 

We subsequently apply the theory to columnar meso- 
phases of charged discotic colloidal particles, a common 
building block for many natural [21] and synthetic clays 
P^24| . The dominant particle interactions stem from 
electrostatic repulsions and the stabilization of columnar 
and other liquid crystal structures can be rationalized on 
entropic grounds alone [25l |26] . Particle density, rather 
than temperature therefore constitutes the chief thermo- 
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FIG. 1: A columnar fluid of stacked parallel soft disks (in 
red) confined in a cylindrical parabolic trap with potential 
Vcir). The backbone of the stacks may adopt a linear (left) or 
buckled (right) configuration. The typical interaction range 
between adjacent disks is quantified by a (Debye) screening 
length . 

dynamic parameter in these lyotropic systems. 

In order to establish a link between the one¬ 
dimensional fluid and the three-dimensional columnar 
meso-phase, we construct a self-consistent mean-field the¬ 
ory by invoking a simple Lennard-Jones cell-theory for 
columnar liquid crystals [?7H^ and identifying the ex¬ 
ternally imposed confining potential acting on the stack 
with an effective internal cell potential imparted by ad¬ 
jacent stacks in the columnar lattice. Expressions for 
the typical amplitude and pitch of the helical undula¬ 
tions as well as the corresponding free energy differences 
between the various buckling modes can then be read¬ 
ily established in relation to the particles density and 
electro-chemical properties of the individual disks via the 
inter-disk potential. We believe that the results of this 
study provide an essential first step in assessing collective 
conformational fluctuations and the possibility of spon¬ 
taneous chiral symmetry breaking in lyotropic columnar 
assemblies. 


II. MODEL 

Let us consider a one-dimensional (ID) fluid of N dis- 
cotic particles at constant temperature T conhned to 
a cylindrical potential trap Vc{r). We focus on dense 
stacks where the average interparticle distance is very 
small so that the disks are assumed to be perfectly 
aligned along the z direction of the laboratory frame 
(see Fig. [^. Their centres-of-mass are constrained to lie 
on a one-dimensional contour parameterized by s(l) = 
{sx{l),Sy{l),Sz{l)} with 0 < I < ^. The contour length 
is given hy ^c = /g and ic > (■ for any contour 

except for a trivial linear one s(Z) = {0,0,/} {^c = 
which corresponds to an ‘unbuckled’ system. Under cer¬ 
tain conditions the linearly stacked column may buckle 
into e.g. a helical shape. 


The central idea of this study is to investigate the opti¬ 
mal buckling shape of the stacked discotic fluid in relation 
to various control parameters such as particle density, the 
choice of the pair interaction, and the external potential. 
The equilibrium shape of the buckled state is not known a 
priori but can be established self-consistently using a sim¬ 
ple variational approach based on classical density func¬ 
tional theory (DFT). To this end we dehne the Helmholtz 
free energy F of a 3D fluid with one-body density p(r) de¬ 
fined as the ensemble average p(r) = (/5(r)) of the micro¬ 
scopic density p{r) = ~ with normalization 

N. The Helmholtz free energy consists of an ideal, exter¬ 
nal held and excess contribution, and is formally given 
by: 

F[p\ = ksT J drp{r) (In A^p(r) - l) 

-f J drp{r)Vcir)+ ^ J drp{r) J dr'p{r')^{r - r'), 

( 1 ) 

with ksT the thermal energy (ks denotes Boltzmann’s 
constant) and A the de Broglie length of each particle. 
While the flrst two contributions are exact, the excess 
term is not. The disk-disk correlations are approximated 
by a simple mean-field term featuring the pair potential 
$(r) between parallel disks. This potential has an intri¬ 
cate directional dependence through the centre-of-mass 
distance vector r which we will specify later on. Let us 
assume the disk centre-of-masses to be distributed uni¬ 
formly along the strictly one-dimensional contour s(/) so 
that we may parameterize p{r) = d/|j9;s||^(r — s(/)) 

in terms of a contour-dependent line density p = Njlc- 
Substituting the parameterization into the free energy 
yields: 

F = -TSid + p[ d/|15zs||W(s(0) 

Jo 

+ ^iyi\\d,s\\l^ d/'||a.s||$(s(0-s(O). ( 2 ) 

Let us consider a helical parameterization s{l) = 
{R cos ql, Rsin ql,l} in terms of radius R and pitch q. 
The contour length of a helix is ic = ^\/l -k {qRY- Since 
a helix has constant curvature, |i9is|| is independent of 
I and the deformation free energy per particle simplifies 
to: 

(3) 

now featuring a linear density p = Nft independent of 
the underlying contour shape. We emphasize that the 
confining external potential is symmetric in the cylin¬ 
drical reference frame of the helix and only depends on 
the lateral extent of the helix via the radius R. Due 
to the one-dimensional constraint imposed on the par¬ 
ticle contour, the ideal entropy Sid cannot immediately 
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be reproduced from its 3D phase-space definition given 
by the first contribution in Eq. Q. In fact, substitu¬ 
tion of the delta-distribution for the one-body density 
into Eq. 0 would lead to a divergence of the free en¬ 
ergy. Nevertheless, a plausible expression for the ideal 
entropy compatible with 3D translational motion in real 
systems can be recovered by factorizing the translational 
degrees of freedom of each disk into a one-dimensional 
contribution representing particle motion along the he¬ 
lical contour and a bi-dimensional one originating from 
infinitesimally small excursions the disk centres-of-mass 
are allowed to make across the plane normal to the con¬ 
tour. This will become apparent later on in this text. We 
will proceed by specifying a suitable form for the inter¬ 
disk potential $. 


III. GENERALIZED YUKAWA POTENTIAL 
FOR CHARGED DISKS 

A natural question that arises is what type of pair po¬ 
tential could give rise to stacked fluids with a stable he¬ 
lical meso-structure? A promising candidate is the ef¬ 
fective electrostatic interaction between charged disks of 
diameter D which in the far-held limit takes the form of 
an effective Yukawa form [3D]: 


$(ri 2 , ui, U2) = uoC( 6 'i; KD)i{62] ^02) ^ 

^12 

with amplitude uq = Z^Xb in terms of the total number 
of elementary charges Z per disk, the Debye screening 
constant k and Xb the Bjerrum length. The angles 9 
featuring in the prefactor ^ denote the disk normal ori¬ 
entation relative to the centre-of-mass distance vector ri 2 
via cos 9 = z- f 12 . The angular prefactor reads in explicit 
form: 


m^D) = 2 


Ii sin 9) 

\kD sin0 


(5) 


IV. FREE ENERGY OF A HELICALLY 
DEFORMED COLUMNAR FLUID 

Let us now attempt to express the directional Yukawa 
potential $ (Eq. 0) in terms of the helical vector s(Z). 
Since the disks normals are all pointing along the z-axis, 
the only relevant coordinate is the centre-of-mass dis¬ 
tance between a pair of particles located at position I and 
I' along the helical contour, i.e., $(Ar) = d)(s(Z) — 
Without loss of generality, the centre-of-mass distance 
Ar = r' — r between two particles located on a helical 
contour can be expressed in the cylindrical coordinate 
frame of one of the particles. Straightforward algebra 
then leads to: 

Ar = R{cos{qAl) — l)f -|- Rsin{qAl)ip + A/z, (6) 

with Al = I' — 1. In order to simplify the orientation- 
dependent prefactor of the Yukawa potential we Taylor 
expand the Bessel function Eq. ([^ for small 9. Some 
manipulation then leads to: 

|Ar| 

where the amplitude depends on the angle 9 between 
the centre-of-mass distance and the disk orientation unit 
vectors aligned along helix axis: 

uo{cos9) = uo (^1 + (l-cosg)^ , (8) 

Up to leading order in q the dot product cos0 between 
the centre-of-mass distance and disk orientations can be 
written as: 


Ar e 

in terms of the contour length e = -|- {qRY- 

amplitude then becomes: 


(9) 

The 


Uo{e) = Uq 
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( 10 ) 


with Ii{x) a modified Bessel function. The Yukawa po¬ 
tential for discotic colloids favors stacked pair configu¬ 
rations over co-planar ones m- As a consequence, the 
inherently electrostatic anisometry drives self-assembly 
into stacked meso-phases, including columnar structures 
m- Conformational changes from a linear to a buckled 
state can be anticipated from a trade-off between the ex¬ 
ternal and excess free energy. While the external free en¬ 
ergy associated with a confining trap potential increases 
upon buckling, the cost is offset by a simultaneous re¬ 
duction of the excess free energy due to the fact that the 
centres-of-mass between adjacent disks will, on average, 
be further apart in an undulated stack. This we will 
explore in detail in the next Section. 


The effective Yukawa amplitude increases with q and R 
and thus clearly disfavors buckling. With this the defor¬ 
mation free energy takes the following form: 


F 

N 


-TSiQ 

N 


+ VciR) + 2“ o (£) 


dAl 


exp[—KAr] 
Ar 


( 11 ) 


with the inter-disk distance Ar depending implicitly on 
the contour parameters A^ via Eq. ([^. The last term 
confronts us with a double contour integral that does not 
converge for unbounded potentials since the disk centre- 
of-masses are forced to reside on a one-dimensional con¬ 
tour. The divergence can be lifted by allowing the disks 
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to perform small spatial fluctuations 5rj_ in the {f, <^}- 
plane perpendicular to the main helical direction z. The 
centre-of-mass distance then takes the form: 


Ar = 6 r+ R{cos{qAl) — l)r + Rsm{qAl)ip + Alz, (12) 


and the contour integral can then be converted to a three- 
dimensional via a straightforward weighted average: 


(13) 

The transverse fluctuations obey a 2D-Gaussian distri¬ 
bution with variational parameter a: 

/G(<5ri)= (14) 


In principle, these fluctuations will also have an effect on 
the angular contribution cos 9. However, their impact is 
expected to be less critical than on the centrally symmet¬ 
ric part and we will retain the form Eq. (10) for simplicity. 


It is expedient to recast the expression in Fourier space in 
terms of the 2D wavevector k. This allows us to express 
the integral in factorized form: 


|uo(e) 


d6r I 


dk 

(2Try 


fG(^r±)e 


— ik.-Sr± 


U(k)P(k), 


(15) 


where H(k) = 47r/(A:^ + k^) is the FT of the spherically 
symmetric Yukawa potential and R(k) a shape factor ac¬ 
counting for the helical shape of the buckled fluid. The 
Fourier transform only acts on the Gaussian displacement 
distribution so that: 

|wo(e) J j^^fG{k)U{k)P{k), (16) 


from the one-dimensional degrees of translational free¬ 
dom along the contour. This leads to a contribution 
—Nksik^Ap — 1) so that the total entropic contribution 
to the free energy takes the following approximate form: 

S'id ~ —.^fcB(ln(pQ(A^/7r) — 3). (19) 

The relevant density variable p = pj^Jl -I- {qRY reflects 
the notion that the available one-dimensional volume 
should increase upon buckling. In the ideal gas limit 
(<& I 0, 14 i 0) a linearly stacked system will triv¬ 
ially maximize its configurational entropy by reducing 
its transverse localization (a decreases) and maximizing 
its fictitious contour length (although the shape need no 
longer be a simple helix). 

Eq. ( |17[ ) cannot be worked out in analytical form for 
arbitrary R. Instead we will resort to a series expan¬ 
sion up to quadratic order in R which should be of suf¬ 
ficient accuracy for weakly buckled fluids. The details of 
this calculation are shown in the Appendix. The generic 
expression for the deformation free energy per particle 
reads: 


^ - \n{pA?a/-K) - 3 - \ {qRf + Vc(i?) 

+ ^uo{e)U{R,q), ( 20 ) 

For notational brevity we shall implicitly use the thermal 
energy kgT as energy unit and k~^ as length unit. An 
approximate expression for the interaction free energy U 
of a weakly buckled helical fluid in the strong localization 
limit (a 1) is derived in the Appendix and reads: 

/ dk - 

—fa{k)U{k)P{k) 

~ (^a)i/2 (1 + R^ap2{q() + 0{{KRf)) , (21) 


where the FT of a Gaussian is itself a Gaussian, namely 
fG{kj_) = The shape factor is formally given by 

a FT of the helical backbone parameterized in Eq. (12): 


.1/1 

P(k) = / (iAZexp[ik • (Ar — i5rj_)], (17) 

J-i/2 

and depends implicitly on the helix shape via Ar. 

Reverting to the ideal entropy in Eq. 0, as of yet un¬ 
specified, we may now associate the transversal Gaussian 
fluctuations with an entropic term of the following form: 

S-^\a) -Nks J dSr±fGiSr±)[\nA'^fG{Sr±)-l] 

~ — AA:B(ln(aA^/7r) — 2), (18) 

indicating a reduction of ideal entropy upon increased 
localization along the contour. As anticipated, S'id di¬ 
verges if the particle centres-of-mass follow a strictly ID 
contour in the limit a —>■ oo. The two-dimensional contri¬ 
bution needs to be supplemented with the entropy arising 


up to quadratic order in R. At this point it is instruc¬ 
tive to broaden our scope a little by making a distinc¬ 
tion between lyotropic systems where Yukawa interac¬ 
tions are repulsive, and thermotropic ones where an at¬ 
tractive Yukawa potential could be used as a rough ap¬ 
proximation for the (7r-)stacking or hydrogen-bonding 
forces between discotic molecules [32l |33] . In the lyotopic 
case pitg > 0 the density is the main thermodynamic 
variable, whereas for thermotropics (pug < 0) the poten¬ 
tial amplitude plays the r ole of an effe ctive temperature 
Teff = |pitg|“^. From Eq. (20) and Eq. (21) we infer that 
an unbuckled fluid of repulsive disks can reduce its excess 
free energy by decreasing a and delocalizing from the 
contour, while attractive disks tends to favor increased 
localization along the contour with F growing more neg¬ 
ative with a. As a consequence, the buckling scenario 
for thermotropics will be quite different from that of ly- 
otropics as we will see shortly. 

Contrary to the ideal contribution, which increases 
with a due to a loss of configurational entropy, the con¬ 
finement free energy must be inversely proportional with 
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particle localization. The free energy penalty associ¬ 
ated with disks moving away from the helical backbone 
can be made explicit by considering a parabolic form 
Vc{r) = vor"^ with vq > 0. The transverse delocaliza¬ 
tion enters via a simple Gaussian average: 

r°° /-Zir 

Vc{R,a) = vo dipfG{Sr±){R +Sr± cos 

Jo ^ Jo 

=»” +i) ^ 


and is inversely proportional to a. In the weak buckling 
limit (i? ^ 1) we assume the equilibrium value of a to 
be unaffected by the undulation so that = 0. 

For the lyotropic case puQ > 0 an approximate analytical 
solution of the minimization condition can be found in 
the limit of strong confinement uq ;S> 1 in which case the 
equilibrium localization strength is governed by a trade¬ 
off between the external and excess free energy. The 
extremum value reads: 

) puo>0, (23) 


showing that the disks become progressively delocalized 
with increasing particle concentration. For the ther¬ 
motropic case the linear fluid is naturally stabilized by 
the strongly attractive face-to-face interactions and the 
confining potential will only have a marginal effect on the 
localization strength. Irrespective of vq > 0, the minimal 
free energy is obtained in the limit a —> oo. In order to 
systematically assess the buckling free energy in powers 
of R we expand the curvature-dependent Yukawa ampli¬ 
tude Eq. (10) up to quartic order in gR 1: 


uo(e) 

Uo 


1 + cf\qRf + c\"'{qRY + 0{R^), (24) 




with = ZlVl6 and cf'^ = (D^ - 48£|2)/1024. It 
is obvious that the orientation-dependent Yukawa am¬ 
plitude in Eq. only makes sense in the context of 
charge-stabilized colloidal interactions where the Debye 
screening length is clearly defined. For the thermotropic 
case the prefactors c„ must have a different origin as 
their amplitude and sign will depend on whether face-to- 
face configurations are energetically favorable compared 
to edge-to-edge ones or vice versa. 

Going back to the lyotropic case we may now formulate 
the free energy difference between the buckled and the 
unbuckled columnar fluid of up to quartic order in the 
buckling amplitude: 


AE(g) 

iVi?2 


1 2 I 1 n2 4 I I / n1/2 

+ 4-^9 +V 0 + —{-KOl) ' 

X + ap2{qf) + c^i^R^q'^ + a^R^Piiqt}^ ■ 

(25) 


In the limit of strong confinement a ^ 1 the free energy 
is fully dominated by the external and excess free energy. 



qH 

FIG. 2: Renormalized buckling radius and associated defor¬ 
mation free energy (on a log scale) for a linear stack of soft 
discotic particles as a function of the optimal pitch values q*£ . 
The energetically most favourable mode is q*i ~ 9.132, with 
corresponding radius R* « 0.333(puo/'Uo)^^^. Higher order 
pitch ‘harmonics’ correspond to increasingly weaker buckling 
radii R* and free energy reductions. 


Minimization with respect to R yields the equilibrium 
value for the buckling radius R*: 


R* CC fiqi) 


puo \ 
Vq ) 


1/3 


< 1 , 


(26) 


in terms of the pitch-dependent function f{q£) = 

+ P 2 {q£))/p 4 {qi) (see Eq. ( |A6| in the Appendix). 
The q values that minimize the free energy turn out to 
be independent of the parameter combination puq/vq and 
follow from f'{q*£) = 0. The corresponding buckling free 
energy is uniformly negative: 


Af(g*) 0+P2{q*£W 

N (/OUo)2/3 Piiq*£) 


(27) 


indicating that helical buckling always leads to a re¬ 
duction of the free energy. On the contrary, for ther¬ 
motropic systems buckling is generally disfavored. Up 
to lowest order in radius it follows that /S.F{q*)/NR? oc 
o?^'^\p 2 {q*£)\/‘FTes 0. This outcome is, of course, not 

surprising given that the attractive stacking interactions 
prevent any deformation of the linear stack at low tem¬ 
peratures. 
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V. CONNECTION TO LYOTROPIC 
COLUMNAR ASSEMBLIES 

As an alternative to the parabolic trap we consider a 
confining potential of the Yukawa form: 


Vyix) = Uo 


1 g-K{cr-x) 

2 a — X 


2^ g-K((T+a:) 

2 a X 


kl < cr, 
(28) 

which can be thought of as a simple effective potential 
acting on each particle along the backbone imparted by 
its neighbouring columns at a typical distance a from the 
target particle. The last term provides the correct off¬ 
set so that VV(0) = 0. Contrary to the parabolic form, 
the Yukawa potential thus features a typical interaction 
range a whose value can either be varied as an inde¬ 
pendent parameter (cf. the spatial extent of an optical 
trap), or it can be coupled self-consistently to the ther¬ 
modynamic state of the columnar assembly. Analogous 
to Eq. (221 the potential is modified by the transverse dis¬ 


placements perpendicular to the column direction. For 
small radii an expansion up to quadratic order is appro¬ 
priate so that: 


Vy{R,a)=uy[R^ + ^+o(?^ 


(29) 


We may now simply repeat the steps outlined in the pre¬ 
vious section. The localization parameter a still obeys 


Eq. (23) but with vq replaced by an effective value as¬ 


sociated with the cylindrical Yukawa potential (denoted 
“Y”), namely uy = UQa~^e~'^ {1 + a + /2). In line with 
the previous section we use k~^ as an implicit length unit. 

A connection with a three-dimensional columnar phase 
can be made by coupling the one-dimensional free energy 
discussed thus far to a simple cell theoretical expression 
for a (hexagonal) packing of columns. Lennard-Jones cell 
theory m tells us that, approximately: 


Qcen(-^) 


drv 

A2 


exp 


Vy{x) ' 

2 kBT 


N 


(30) 


Within this framework the line fluids are localized in N 
cylindrical cells centered on the sites of a fully occupied 
lattice of some prescribed symmetry. Each particle expe¬ 
riences a potential energy Vy (R) generated by its nearest 
neighbors. The term between brackets is therefore equiv¬ 
alent to the lateral ‘free volume’ each particle explores 
within its cell. In its simplest version, the theory presup¬ 
poses each cell to contain only one stacked fluid behaving 
independently from its neighbors pTIEH] . Moreover, the 
free volume term does not reflect the precise nature of 
the 2D Bravais lattice of the columnar phase which we 
assume to be a simple hexagonal arrangement. Using 
Eq. (29) and some algebra then yields for the partition 


sum: 


N 


Qce\liN) 


I — e 


Uyh? 


(31) 


Single occupancy of each cell implies that the three- 
dimensional density g = N/V is linked to the linear 
density via g = Recalling that F = —fc^TlnQ 

we obtain the free energy per particle of the colum¬ 
nar phase by superimposing the contributions from the 
stacked fluid and the cellular confinement (ignoring irrel¬ 
evant constants) P5][55]: 


f{q) 


\n{ga‘^a) + + ^uo(7ra)^/^ -I- Inuy 

iV la 1 

- In f 1 - e 


1^1. 


(32) 


In the limit of very strong confinement we expect a 
so that Uy ^ uo/a^ 3> 1. Minimization of the colum¬ 
nar free energy with respect to a yields, as expected, an 
inverse proportionality of the localization strength with 
density and amplitude, i.e. a ~ 

Subsequent minimization of the free energy with re¬ 
spect to cr gives a similar proportionality but with dif¬ 
ferent exponents, namely a ^ crg/upp^ with constant 
CTo = 70304/277r. With this we can establish an explicit 
relation between the density g of the columnar phase and 
the localization strength, namely: 


aoul°g^, 


(33) 


with prefactor ag = (4/7r)^/^crg Similarly, the effec- 
tive amplitude of the self-consistent cell potential reads: 


. 10^6 


Uy 


(34) 


In order to get a feeling of the magnitude of the variables 
above we may identify uq « Z'^kXb with Z ^ 1 the 
number of bare or effective elementary charges per disk 
[301 [M] and Xb the Bjerrum length (which is typically 
about I nm in aqueous conditions at room temperature). 
The effective columnar packing fraction can be estimated 
from (j) ~ N/V so that g ^ /D)'^ < 1. The 

typical buckling radius then scales as: 








(35) 


and the deformation free energy associated with the pitch 
optima (see Fig. scales as follows: 


AE(g*) 

N 


,20 12 
ig g 


{l+P 2iq*m^ 


< 0 . 


(36) 


For infinitesimally weak buckling deformations, a Taylor 
expansion of Eq. (25) up to the lowest order contribu¬ 


tion in the helix parameters gives an expression which 
could be interpreted as the typical bend or twist elastic 
modulus of a single column: 


AF 

N{qRy 


,10 




24^3 


(37) 
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The modulus is negative because internal stresses along 
the backbone lead to buckling instabilities and ‘floppy’ 
modes whose amplitude increases quadratically with 
length £. We finally wish to point out that the buck¬ 
ling radius and free energy depend implicitly on the De¬ 
bye screening length. Rearranging the scaling expressions 


revea 


Eqs. (35) and (36 1 in terms of the (fixed) disk diameter D 
, s that R*YD cx {D/k~^Y AF/N oc 


and shows that the buckling characteristics depend sen¬ 
sitively on the range of the stacking potential. 


VI. CONCLUDING REMARKS 

We have proposed a microscopic theory to probe heli¬ 
cal buckling deformations of a columnar fluid of discotic 
particles. Contrary to the more commonly studied case 
of buckling under an external load, the stresses that sta¬ 
bilize helical confirmations of the columnar backbone are 
purely internal and are transmitted through a strongly 
repulsive stacking potential between adjacent disks lo¬ 
cated along the backbone. The approach is particularly 
suited to investigate conformational fluctuations in ly¬ 
otropic columnar assemblies of e.g. clay platelets at low 
salt concentrations where the effective charge-mediated 
particle interactions are strongly repulsive. Using a sim¬ 
ple mean-field theory based on microscopic principles we 
derive scaling relations for the optimal pitch, which sets 
the typical buckling length scale, and amplitude in rela¬ 
tion to the density of the columnar mesophase and the 
strength of the electrostatic interactions. Our findings 
indicate that helical deformations of stacked fluids can 
be stabilized by long-ranged soft interactions. 

Our analysis is based on the assumption that the 
columnar deformations are independent from each other. 
However, we consider a detailed knowledge of the single¬ 
column buckling behaviour as an essential precursor to 
tackling the more complicated problem of collective con¬ 
formational fluctuations. In this context, we plan to ex¬ 
tend the present model such as to include chiral inter- 
columnar couplings, e.g., by introducing a weakly chiral 
self-consistent cell potential, in an effort to understand if 
and how helical conformational fluctuations of an individ¬ 
ual column proliferate to larger length scales. This could 
open a route to studying spontaneous chiral symmetry 
breaking and chiral amplification in columnar assemblies 
of achiral discotic particles under the influence of inter¬ 
nal stresses generated by repulsive long-ranged stacking 
potentials. The results of the current model are expected 
to facilitate the construction of coarse-grained models 
for studying the chiral intra-columnar meso-structure 
of columnar liquid crystals. These models could pro¬ 
vide a useful intermediate route between fully particle- 
resolved simulations |35L I36j and phenomenological con¬ 
tinuum theories that have been fruitfully employed to 
study columnar bundles of helical fibres m and defects 
in chiral columnar phases [3S]. Research efforts along 
these lines are currently pursued. 


Appendix A: Excess free energy of helically stacked 
soft disks 

Let us decompose the wavevector in a 2D transver¬ 
sal and longitudinal component so that / dk = 
f dtp f dkj_k± f dk\\ and -I- fcjj. The Fourier inte¬ 

gral Eq. can then be expressed as follows: 

/•OO 11 p D7 

ld{R,q) = -^kjcik^) J -^U{k)P^{k), (Al) 

in terms of the the angle-averaged shape factor 

J-ljl 

(A2) 

The above integral can be greatly simplified by pre¬ 
averaging the pair potential over the longitudinal wave- 
vectors, an approximation that does not lead to qualita¬ 
tive error |40] : 



U{R,q) « ^fcu/G(fcj.)aA(fcu)(U(fc))||. (A3) 

Taking the FT of the Yukawa potential U(k) = 47r/(fc^ -|- 
K^) gives: 

/ OO I - 

^U{k) = 27r/^kl+K^, (A4) 

which now depends only on the transverse Fourier com¬ 
ponent. In the weak buckling limit (i? 1) we may 

expand up to quartic order in the buckling amplitude 
R. Retaining only the real-valued contributions yields: 

Ptp{k±) = 1-1- - —^d-p 2 (gf) -I- ^ 32 ^ PiidP) + 0{R^), 

(A5) 


in terms of the following pitch-dependent functions: 
P2iq£) = joiqi/^) - 1 

Piiqi) = 3 - Ajo{qi/2) + jo{qi), (A 6 ) 

with jo{x) = s\Tix/x a spherical Bessel function. Note 
that the helix form factor is an even function of q and 
does not discriminate between a right-handed (g > 0 ) or 
left-handed {q < 0) helix sense. The integrations over the 
transverse wavenumber k± can be performed analytically. 
In the limit of asymptotically large a it sufhces to retain 
only the leading order powers in a: 

r°° Hk, ^ 

X ~C 2 „(a), (A7) 

with Co = (ttq:)^/^, = 27r^/^Qf^/^ and C 4 = . 

Combining terms we obtain the following expression for 







the interactions free energy in the strong localization 
limit for arbitrary order in the pitch q: 
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